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ABSTRACT 



Deeply embedded and at distances of several kiloparsecs, massive young stellar objects (MYSOs) present numerous challenges for 
observation and study In this work, we present spatially-resolved observations of one MYSO, AFGL 4176, together with survey and 
literature data, ranging from interferometric observations with VLTI/MIDI in the mid-infrared, to single-dish Herschel measurements 
in the far-infrared, and sub-millimeter data from APEX. We consider this spatially-resolved, multi- wavelength data set in terms of 
both radiative transfer and geometric models. We find that the observations are well described by one-dimensional models overall, 
but there are also substantial deviations from spherical symmetry at scales of tens to hundreds of astronomical units, which are 
revealed by the mid-infrared interferometric measurements. We use a multiple-component, geometric modeling approach to explain 
the mid-infrared emission on scales of tens to hundreds of astronomical units, and find the MIDI measurements are well described by 
a model consisting of a one-dimensional Gaussian halo and an inclined (6 = 60°) circumstellar disk extending out to several hundred 
astronomical units along a position angle of 160°. Finally, we compare our results both with previous models of this source, and with 
those of other MYSOs, and discuss the present situation with mid-infrared interferometric observations of massive stars. 

Key words, stars: formation - techniques: interferometric - techniques: high angular resolution - radiative transfer - stars: individual: 
AFGL 4176 



1. Introduction 

The class of objects known as massive young stellar objects 
(MYSOs) are thought to represent massive stars at the earliest 
phases of evolution. The inherently large distances (several kpc) 
and high extinctions (A v ^ 20 mag) to typical MYSOs compli- 
cate attempts to c haracterize these obje cts observationally, and 
many studies (e.g. lGrave & Kumarf20 09) rely on fitting the spec- 
tral energy distribution (SED) to untangle the underlying physics 
involved in the formation of massive stars. 

However, it has been shown that pure SED fitting, 
i.e. without making use of any spatia l information is a 
highly degenerate approach, at best (e.g. iThamm et al.l Il994t 
iMen'shchikov & Henningi[T997h . Different physical scales and 
temperature regimes of the circumstellar material can be re- 
solved using a variety of techniques at different wavelengths, 
ranging from weakly-resolved single-telescope measurements, 
spectroscopic (spectroastrometric) measurements, and interfer- 
ometric measurements, both at near-/mid-infrared and (sub- 
millimeter wavelengths. The temperatures and regions thus 
probed range from ~ 10 3 - 10 4 K gas within the dust sublima- 
tion radius using spectro astrometry of molecular gas lines (e.g. 
IWheelwright et aUboiOl) . ~ 10 3 K dust at the inner dust disk ra- 
dius using near-infrared interferometry (e.g. iKraus et ai]|2010l) . 
~ 10 2 - 10 3 K dust farther out in the disk (and possib ly into out- 
flow cones) using mid-infrared interferometry (e.g. iLinz et all 
120091) . and the broad distribution of cold ~ 10 1 - 10 2 K dus t 
using millimeter- wave interferome try (e.g. iBeuther et al.ll2007l) 
and single-dish measurements (e.g. lBeltran et al.ll2006l) . 



In this paper, we present observations and the results of 
multi-wavelength modeling of one particular massive young 
stellar object, AFGL 4176 (also known as IRAS 13395-6153, 
G308.9+0.1). This object was identified as a candidate MYSO 
based on the similarity of the far-infrared spect rum to that of 
the Becklin -Neugebauer (BN) object in Orion dHenning et al.l 
fT98lfT990h . The class of so-called BN-type objects are lumi- 
nous (^ 10 3 Lq) sources, extremely red (usually with no opti- 
cal counterpart), with a broad peak in emission (flat spectrum in 
vF v ) between 50 and 200 jim and deep absorption due to silicates 
at 10 fim. They are generally interpreted as young (possibly still 
forming) massive stars, surrounded by a thick dust envelope, de- 
spite the un clear natur e of the prototype object Orion BN (see 
the work by lTanl (120041) . for example). 



At radio wavelengths, us ing the Molonglo O bservatory 
Synthesis Telescope (MOST), ICaswell et afl (1 19921) reported a 
small-diameter source (<0f2) roughly coincident with the IR 
source AFGL 4176. The integrated flux densities are 0.31 Jy 
at 843 MHz and 0.40 Jy at 1415 MHz, confirming the emis- 
sion is thermal bremsstrahlung radiation arising from an H II re- 
gion. Using the Aus tralian Telescope Compact Array (ATCA), 
iPhillips et al.l (1 19981) reported compact continuum emission at 
8.6 GHz spanning about 10", with a collecti on of 6.7 GHz 
metha nol masers on the northern edge. Later, lEllingsen et al.l 
(120051) showed the 8.6 GHz continuum emission extends to 
about r. Finally, th e 1.2 mm dust continuum was mapped by 
iBeltran et al.l (120061) . and comprises a roughly circular, ~ 40" 
clump at the same location with an estimated mass of 1 120 Mq 
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at 5.3 kpc, and a fainter clump (estimated mass 177 Mq at 
5.3 kpc) about V.5 to the east. 

In the near infrared, AFGL 4176 is a point source for seeing- 
limited obse r vation s. Using speckle interferometry, however, 
iLeinert etaD (1200 ll) found a faint extended halo in the K band, 
with a size of Y.'6 ± 074. The positions from the Two Micron 
Al l Sky Survey (2MA SS) and the methanol maser s repo rted 
by iPhillips et al] d 19981) are within a/7. iDe Buizeri (l2003h de- 
tected a large number of shocked H2 emission blobs around 
AFGL 4176 by means of narrow-band NIR imaging. However, 
their spatial distribution does not reveal a preferred direction for 
a potential jet, and in general, the association of AFGL 4176 
with a jet or outflow remains unclear. The source was also ob- 
served in the mid- through far- infrared by the space missions 
IRAS, MSX, ISO, Spitzer and Herschel. We note that the coordi- 
nates listed in the IRAS Point Source Catalog appear to be erro- 
neous when compared with the MSX 12 fim and Spitzer 24 jim 
images, and adopt the 2MASS position of a = 13 h 43 m 0P.70 
6 = -62°08 , 51 , /23 (J2000) as the location of the infrared source. 
The ISO spectrum shows deep silicate absorption features at 9.7 
and 1 8 yum, as well as absorption due to water ice at 3 fim, CO2 
ice at 4.3 yum, and emission lines of [Nell] at 12.8 yum; [S III] at 
12.0, 18.7 and 33.5 yum; [01] at 145.5 yum and possibly 63.2 //m; 
and [C II] at 157.7 /im. 

The distance to AFGL 4176 r emains poorly constr ained. 
Using both the velocity curve of iBrand & Blitz! (1 1993b with 
updated par ameters from ILevine et al.1 (l2008l) .~and the veloc- 
ity curve of iReid et al.1 d2009h , we find that the combination of 

the source's LSR velocity of 51 km s" 1 (e.g. iFontani et al.l 

De Buizer et ID 12009) and galactic coordinates are for- 
biddenj. Consequently, a kinematic distance cannot be deter- 
mined for AFGL 4176 using these rotation curves. Despite 
this, several authors (e. g . iBeltran et al .120061: IFontani et al.ll2 005 ; 
iGrave & Kumarl 120091) r eport kinematic distances of around 
5 kpc. ISaito et all (l200ll) showed that the spread in velocities 
of C 18 molecular clumps in the vicinity is high (~ 30 km s" 1 ), 
meaning that kinematic estimates of the distance to this object 
should be highly uncertain, even after improvement of the galac- 
tic structure parameters. To date, other means of the distance 
measurements (e.g. maser parallax measurements) have not been 
applied to AFGL 4176. Because of this uncertainty, we consider 
the observational data in terms of both "near" (3.5 kpc) and "far" 
(5.3 kpc) distances, corres ponding to the edg es of the large C 18 
complex as determined by lSaito et al.l (l200ll) . 

The present paper is similar in scop e to previous in - 
frared interf e rometr i c studies of MYSOs b y iLinz et al J (l2009h : 
iFollert et al] d2010h : Ide Wit et al.1 d2010h : iKraus et al.1 0201 Oh . 
However, in this work on AFGL 4176, we go beyond tradi- 
tional SED modeling techniques and use a far more compre- 
hensive approach, including not only flux values, but also spa- 
tial information at mid-infrared, far-infrared, and sub-millimeter 
wavelengths directly into the fitting process. Such an approach 
helps to break the degeneracies inherent in modeling the SEDs 
of deeply-embedded sources, and provides structural informa- 
tion on scales ranging from tens of AU up to a parsec. We begin 
by presenting and discussing new observations of AFGL 4176, 
then go to one-dimensional radiative transfer modeling of the 
entire data set. Next, we use a geometric component analysis of 
the mid-infrared interferometric data to explore deviations from 



spherical symmetry, which we interpret as strong evidence for 
a circumstellar disk around this object. Finally, we compare the 
results of our study with those of previous studies of this object, 
as well as studies of other MYSOs. 



2. Observations and data reduction 

2.1. Mid-infrared interferometry 

AFGL 4176 was observed with the two-telescope mid-infrared 
interferometer MIDI on the Very Large Telescope Interferometer 
(VLTI) of the European Southern Observatory on Cerro Paranal 
in Chile. Observations were conducted in 2005-2007 and 2011- 
2012 with both the 8.2 m Unit Telescopes (UTs) and the 1.8 m 
Auxiliary Telescopes (ATs) as part of guaranteed-time obser- 
vations (proposal IDs 074.C-0389, 075.C-0755, 076.C-0757, 
077.C-0440, 078.C-0712 and 084.C-1072). In total, 37 points 
in the uv plane were obtained for projected baselines ranging 
from 5 to 62 meters and a variety of position angles, with six 
of the measurements obtained with the UT telescopes, while the 
remaining 31 were obtained with the AT telescopes. This repre- 
sents, by far, the largest collection of long-baseline mid-infrared 
interferometric observations presented for a MYSO to date. A 
summary of the observations is shown in Table [T] and the uv 
coverage is shown in Fig.[T] For clarity, the individual observa- 
tions have been color-coded by projected baseline length: purple 
indicates measurements with the longest baselines (60-64 m); 
red (36-49 m) and green (25-32 m) show those with interme- 
diate baselines; blue (13-16 m) and orange (5-11 m) mark the 
shortest baseline measurements. 

We used the prism as the dispersive element, resulting in 
spectrally-resolved visibility amplitudes and differential phases 
with a resolution of A/AA = 35 in the wavelength range 8 - 
13 fim. All observations were carried out in the HIGHSENS 
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1 We note that the tangential solution of these velocity curves at the 
coordinates of AFGL 4176 occurs for V LS r « -49 km s _1 and corre- 
sponds to about 5.3 kpc, which may explain the kinematic distances 
reported by other authors. 



Fig. 1. uv coverage of MIDI observations. North is up, east is 
left. The positions have been color-coded by projected baseline 
length, and the dashed lines indicate the mean projected base- 
lines for each color (orange - 8.1 m; blue - 15.2 m; green - 
28.5 m; red - 42.5 m; purple - 61.7 m). 
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Date/time 


Telescopes 


Proj. baseline 


Position angle 


Avg. visibility 


(UTC) 




(m) 


(deg.) 




2005-03-02 08:08 


UT3-UT4 


61.6 


112.1 


0.02 


2005-06-24 04:54 


UT3-UT4 


61.8 


169.3 


0.02 


2005-06-25 22:50 


UT1-UT2 


49.3 


19.3 


0.07 


2005-06-26 02:05 


UT1-UT2 


41.9 


49.8 


0.08 


2005-06-26 03:27 


UT1-UT2 


36.3 


62.4 


0.04 


2005-06-27 03:31 


UT2-UT3 


30.3 


85.6 


0.08 


2006-02-22 05:46 


E0-G0 


16.0 


41.9 


0.30 


2006-04-18 08:01 


D0-G0 


26.5 


121.4 


0.06 


2006-04-24 07:30 


E0-G0 


13.4 


119.4 


0.31 


2007-02-11 06:24 


E0-G0 


16.0 


40.5 


0.23 


2007-02-12 06:18 


GO-HO 


31.9 


40.0 


0.13 


2007-03-15 07:08 


GO-HO 


31.2 


77.7 


0.11 


2007-03-16 09:13 


GO-HO 


28.3 


106.8 


0.07 


2007-04-18 04:55 


EO-GO 


15.6 


77.8 


0.20 


2007-05-24 01:51 


EO-GO 


15.8 


68.7 


0.27 


2007-06-19 02:35 


GO-HO 


28.9 


101.0 


0.06 


2007-07-03 01:31 


EO-GO 


14.6 


98.9 


0.28 


2011-05-09 02:15 


B2-C1 


9.8 


21.7 


0.47 


2011-05-09 02:49 


B2-C1 


9.6 


27.0 


0.46 


2011-05-09 03:14 


B2-C1 


9.5 


31.0 


0.49 


2011-05-09 03:29 


B2-C1 


9.4 


33.4 


0.52 


2011-05-09 04:54 


B2-D0 


25.9 


46.4 


0.18 


2011-05-09 05:27 


B2-D0 


24.7 


51.6 


0.16 


2011-05-09 08:28 


B2-C1 


5.3 


80.9 


0.61 


2011-05-10 02:23 


B2-D0 


29.3 


23.5 


0.20 


2011-05-10 03:46 


B2-D0 


27.7 


36.7 


0.19 


2012-02-01 09:04 


A1-B2 


10.7 


103.5 


0.51 


2012-02-24 05:20 


DO-HO 


63.8 


37.3 


0.06 


2012-02-24 05:30 


DO-HO 


63.8 


39.8 


0.06 


2012-02-24 09:29 


DO-HO 


59.9 


92.0 


0.05 


2012-03-26 09:05 


Il-KO 


32.4 


38.6 


0.11 


2012-03-26 09:41 


Il-KO 


30.8 


43.1 


0.10 


2012-03-27 05:21 


Il-KO 


38.0 


9.1 


0.10 


2012-03-27 05:26 


Il-KO 


37.9 


9.7 


0.09 


2012-03-27 05:34 


Il-KO 


37.8 


10.9 


0.09 


2012-03-27 07:09 


Il-KO 


36.2 


23.8 


0.10 


2012-03-29 04:08 


Il-KO 


38.2 


179.9 


0.09 



Notes. For clarity, observations are color-coded by projected baseline length, using the same color scheme as Fig.[T] 



mode, meaning the correlated flux and the photometric measure- 
ments are done separately. For a d etailed description of this ob- 
serving procedure, see the work by lLeinert et al.l (120041) . In gen- 
eral, every science measurement was preceded and/or followed 
by a calibration measurement of a bright single star of known 
diameter. 

The raw data were reduced using the MIA+EWS package 
version 1.7.1. Masks were generated for each interferometric 
measurement from the corresponding photometric data. From 
the calibrator measurements taken throughout the night, we 
computed the mean sensitivity function, relating counts/s to Jy of 
correlated flux. The uncertainty in this sensitivity function was 
taken to be the standard deviation over all the calibrator mea- 
surements in the night, and is typically on the order of 10%. 

Because the quality of the photometric data from the ATs is 
generally poor, we assume AFGL 4176 to be non- variable and 
use an average total spectrum from the UTs for computing the 
visibility amplitude from the correlated flux. Uncertainties in the 
final visibilities were calculated using standard error propagation 
formulae. 



2.2. Far-infrared imaging 

The AFGL 4176 region was covered by far- infrared observation s 
taken by the Herschel Space Observatory dPilbratt et al.l 120 1 pj) 
within the framework of the Hi-GAL survey dMolinari et al.l 
12010b . In this paper, we present bolom eter data at 70 and 160 yt/m, 
obtained with the PACS instrument dPoglitsch et al.l l2010h on 
August 16, 2010. These data are scan maps with a scan rate of 
60"/s, which slightly degrades the spatial resolution compared 
to the theoretical diffraction limit. After onboard preprocessing, 
the data are transmitted with a downlink rate of 5 frames/s for 
the 70 jjm filter, and 10 frames/s for the 160 fim filter. The final 
resulting spatial resolution for the 70 jim data is about 9 / /2, and 
1276 for the 160 jim data, although the true PSF is asymmet- 
ric. For a more detailed description of the observing procedure 
used in the Hi-GAL program, we refer to t he description of the 
observing strategy by lMolinari et al.l (120101) . 

We extracted the level-0 data from the Herschel Science 
Archive of the relevant 2° x 2° strip covering AFGL 4176 
(OBSIDs 1342203085 and 1342203086). Despite the occurrence 
of very bright far-infrared sources in the Hi-GAL fields, these 
observations were performed with the normal high-gain settings 
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for the PACS bolom eters The level-0 data were processed with 
the HIPE software dOttll2010h version 8.0.3178 up to the so- 
called level- 1, where internal data units have been translated to 
Janskys, and outliers (i.e. glitches) in the data and saturated pix- 
els have been removed or masked. 

In addition to these standard procedures, we applied a task 
for eliminating the electronic cross-talk within the bolometers, 
which can cause the appearance of ghosts close to strong com- 
pact sources. Furthermore, we use a module recently developed 
by the instrument team which applies a non-linearity correc- 
tion to the data. This is especially important for our purpose, 
since in the high-gain mode used for the Hi-GAL program, very 
bright sources will trigger a non-linear response in the bolometer 
behavior, resulting in an altered PSF and decreasing measured 
fluxes by 10-15%. 

The removal of l/f noise and effects resulting from tem- 
perature drifts within the bolometers, as well as the final map- 
ping, was done usin g the version 7 of the Scanamorphos pro- 
gram (lRousselll2012l submitted to PASP), for which we used the 
"galactic," "noglitch" and "parallel" options. Final maps were 
produced with a pixel scale of 3" 2 for both filters, which is the 
native plate scale for the 70 fim bolometers, and half the native 
plate scale for the 160 fim bolometers. 

2.3. 870 jum imaging 

AFGL 4176 was also imaged at 870 jim as part of the APEX 
Telescope Large Area Surv ey of the Galaxy (ATL AS GAL). We 
refer to the publication by ISchuller et al.l (l2QQ9h regarding the 
general outline of this sub-millimeter continuum survey, and the 
technical details of the observations and the data reduction. In 
the present work, we use a map from the re-reduction of the 
data in 2011, done by the F. Schuller, using refined methods 
for thermal drift removal within the bolometer reduction soft- 
ware BOA (ISchuller et al.ll2009l) . Furthermore, the treatment of 
strong sources within the iterative process of the data reduction 
has been improved. Finally, the data have been mapped onto a 6- 
arcsec grid, u sing version 2.18 of the mapping software SWarp 
(lBertinll2(5TQb . which prevents the formation of spurious Moire 
patterns still present in earlier ATLASGAL maps. 
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Fig. 2. Above: Acquisition image of AFGL 4176 from MIDI ob- 
tained on the 8.2 m UT3 telescope using the SiC filter, which 
has a central wavelength of 11.8 /mi and a width of 2.32 jim. 
Below: Comparison of a cut through the acquisition image of 
AFGL 4176, and an analogous image of the unresolved calibra- 
tor starHD 102461. 



3. Results 

We show an AO-corrected acquisition image of AFGL 4176 ob- 
tained with the 8.2 m telescope UT3 on Mar. 2, 2005 in the 
upper panel of Fig. [2 The image was obtained using the SiC 
filter, which has a central wavelength of 1 1.8 jim and a width of 
2.32 fim. In the lower panel, we compare a cut through the acqui- 
sition image of AFGL 4176 with an analogous cut through the 
acquisition image of the unresolved calibrator star HD 102461, 
which was observed immediately after the science source. From 
the similarity of the spatial cut through the image of AFGL 4176 
to that of the calibrator, we conclude that AFGL 4176 is at most 
marginally resolved in diffraction-limited 8 m-class telescope 
images at Af-band wavelengths. 

In Fig. [3j we show the average Af-band spectrum of 
AFGL 4176 from the 8.2 m UT telescopes, used to derive the 
visibility amplitude from the correlated flux measurements. The 
visibility amplitudes as a function of wavelength are shown in 
Fig.|4l where the measurements have been ordered by projected 
baseline length, and are presented using the same color scheme 
as in Fig. [TJ Finally, in Fig. we show the differential phases 
measured with MIDI for the measurements with B > 36 m (the 
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Fig. 3. Average Af-band spectrum measured with MIDI on the 
UTs. 



4 



Paul A. Boley et al.: On the massive young stellar object AFGL 4176 





(J.D 




(J. 4 




0.2 




0.0 




0.O 




0.4 




0.2 




0.0 




6 




4 




0.2 




0.0 




U.D 


(D 


4 






1 1 


0.2 






E 


0.0 




U.o 


>^ 
-<— » 


!q 


U.Z 




0.1 


> 




0.0 




0.O 




U.Z 




0.1 




0.0 




0.O 




O.z 




0.1 




0.0 




0.3 




0.2 




0.1 




0.0 




n 1 r 



B = 5.26,m =,80.9° 



B f= 9.38,m 6> =,33.4° 




4^4 



9.48 m 6> = 31° 





9 10 11 12 9 10 11 12 9 10 11 12 

Wavelength (/im) 



Fig. 4. The complete set of visibility amplitudes measured with MIDI. The plots are sorted by projected baseline length, indicated 
as B in the inset text, together with the position angle 6. The color coding is the same as for Fig.Q] 
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Fig. 4. Continued. 



differential phases at shorter baselines are not shown, as they are 
essentially zero). 

We note the clear presence of a zero crossing in both the 
visibility amplitude and differential phase for the two measure- 
ments with 5 = 61.6 and 61.8 m, indicated with a vertical dashed 
line in Figs. |4] and [5] We previously interpreted this as evidence 
for ci rcular symmetry at the smallest scales probed (iBoley et al.1 
l201ll) . In order to test this hypothesis, we targeted this spatial 
frequency for follow-up observations with the ATs at other po- 
sition angles on Feb. 24, 2012 (the remaining "purple" measure- 
ments in Figs. [T]|4] and 0. However, while a non-zero differen- 
tial phase is measured, a zero crossing is clearly not present in 
these follow-up observations, implying deviations from spheri- 
cal symmetry at the smallest angular scales probed by MIDI. 

In Fig. [6j we show the images at 70 and 160 fim (Hi-GAL) 
and 870 /im (ATLASGAL). The approximate beam size for each 
image is 97 '2, 1276 and 1972, respectively, and is shown as a 
circle in the lower left-hand corner of each image. The position 
of the infrared source is indicated as a cross in each image. The 
position of the emission peak in the Hi-GAL data corresponds 



to that of the 2MASS source to within the pointing accuracy of 
the observations. In the ATLASGAL data, the emission peak lies 
~ 4" south of the IR position, although this is also within the rms 
pointing accuracy of the ATLASGAL observations. 

At the Hi-GAL wavelengths, within a 36" aperture, we mea- 
sure flux densities of 2382 Jy and 1383 Jy at 70 yum and 160 yum, 
respectively. While the absolute flux calibration on these data 
is accurate to within 5% (Herschel PACS Report PICC-ME- 
TN-037, Apr. 12, 2011), there is significant, non-uniform back- 
ground emission in the region which is difficult to account for. 
Therefore, these values should be treated as upper limits, as it is 
not clear where the background begins to dominate and the ex- 
tended emission associated with the source stops. The 870 jim 
emission is also quite extended; within a 19" aperture, which 
covers just the central emission peak, we measure an integrated 
flux density of 2.86 Jy, while in a larger 72" aperture, covering 
the extended emission, we measure an integrated flux density of 
10.6 Jy. 

In Table [2] we list the photometric data and effective aper- 
ture sizes from previous works and surveys which we consider 
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Fig. 5. Differential phases measured with MIDI for uv points with projected baselines longer than 36 m. The vertical dashed line for 
the two measurements with B = 61.6 and 61.8 m shows the location of the phase flip, indicative of a zero crossing in the visibility 
amplitude. 



in modeling our source. Unless noted otherwise, we compare 
our model data (see Sec. 14. lb to the observations by considering 
the flux contribution within a Gaussian aperture with the same 
FWHM as the observational data. 

In Fig. [3 we show H- and A^-band images of AFGL 4176 
and its immediate surroundings from the first data release of the 
ESQ public surv ey VISTA Variables in the Via Lactea (VVV; 
ISaito et all 120121) . These images were obtained on the 4 m 
VISTA telescope on Cerro Paranal, and surpass 2MASS both 
in spatial resolution and sensitivity. The seeing for the /, H and 
K s images presented here is 0782, 0786 and (J/94, respectively. 
AFGL 4176 is present as a bright, compact source, and saturated 
in both the H and K s images. 

Finally, the region around AFGL 4176 was mapped in sev- 
eral molecular lines as part of the Millimeter Astronom y Legacy 
Team 90 GHz (MALT90) survey dFoster et alJl2011b with the 
Mopra 22 m telescope. Of the 16 spectral lines which were 
mapped as part of the survey, HCO+ was the brightest in 



AFGL 4176. The HCO+ line data, shown in Fig. H have a 
spectral resolution of ~ 0.11 km s" 1 . This line is a tracer of 
dens e molecular gas, an d has been used as an outflow indicator 
(e.g. lHofner et al.ll200ll) . However, in the case of AFGL 4176, 
the spectrum is very well approximated by a Gaussian (FWHM 
3.5 km s" 1 ), and shows no signs indicating the presence of an 
outflow. 



4. Analysis 

4.1. One-dimensional radiative transfer modeling 

As shown in Sec. [3j AFGL 4176 shows evidence for devia- 
tions from spherical symmetry only at the scales probed by our 
mid-infrared interferometric observations, which correspond to 
scales of roughly 20 to 200 mas. Further more, no outflow c om- 
ponent has been revea led in searches by iDe Buiz er (2003) and 
iDe Buizer et al.l ([2009), nor in the HCO+ data presented in this 
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Fig. 6. The 70 and 160 fim images obtained from the Hi-GAL survey with the Herschel Space Observatory, and the 870 fim image 
from the ATLASGAL survey with the APEX telescope. The circle in the lower left corner of each image indicates the beam size. 
The two Hi-GAL images are shown with logarithmic scaling of the intensity, while the ATLASGAL image is shown with linear 
scaling. 
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Fig. 7. H- and ^-band images from the VVV survey. AFGL 4176, indicated with a white cross, is saturated in both the H and 
K s images. The circle in the lower left corner of each image indicates the seeing, as measured from the guide star by the autoguider. 
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Fig. 8. Spectrum of HCO+ from the MALT90 survey. The ob- 
servations are shown with the gray line, while a Gaussian fit is 
shown with the black line. 



work (Fig. [8]). Since there is no evidence for asymmetry or pre- 
ferred geometry on scales larger than a few hundred AU, we be- 
gin with one-dimensional models to explore the overall structure 
of this source. 



We model AFGL 4176 using the one-dimensional radiative 
transfer code MoDust, which was developed by A. de Koter and 
J. Bouwman. MoDust use s the Feautrier m ethod to solve the ra- 
diative transfer eq uation (lFeautrierlll964l). and has be e n use d 
previously by e . g. iMalfait et al.1 (11999b : iKemper et al.1 (l200ll) : 
iBouwman et al.l (120011) . Using MoDust. the dust temperature 
structure is iteratively solved for a given density structure and 
dust properties. Different opacities can be given for different 
temperature regimes, thus making it possible to dynamically de- 
termine the distributions of different types of dust throughout the 
model. 



Using this code, we produce radial intensity profiles of the 
source image on the sky at all wavelengths for which we have 
observations, i.e. from 1.25 yum to 1.2 mm. From these radial 
intensity profiles, we calculate aperture-corrected flux values for 
the SED, PSF-convolved images for comparison with the Hi- 
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Table 2. Photometric data 



Wavelength 


Flux density 


Aperture 


Reference 


(jim) 


(Jy) 


(") 




1.25 


(2.664 + 0.139) x 10" 3 


4 


(1) 


1.65 


0.123 + 0.004 


4 


(1) 


2.20 


1.760 + 0.04 


4 


(1) 


9.0 


173.2 + 1.33 


11 


(2) 


12.0 


250.70 + 12.535 


— 


(3) 


18.0 


434.2 + 13.3 


11 


(2) 


25.0 


601.60 + 30.080 


— 


(3) 


60.0 


2701.00 ±243.09 


— 


(3) 


65.0 


1506.0 ± 434.0 


27 


(2) 


70.0 


2382.2 


36 


(4) 




qaci 1 i ono r\ 
o4^.1 ± zUo.U 


11 


(2) 


100.0 


3566.00 ± 534.90 




(3) 


140.0 


1437.0 ± 127.0 


44 


(2) 


160.0 


1939.0 ± 293.0 


44 


(2) 


160.0 


1383.4 


36 


(4) 


870.0 


2.86 


19 


(4) 


870.0 


10.6 


72 


(4) 


1200.0 


3.86 


39 


(5) 



Re ferences. (1) 2MA SS; 
f5) lBeltraneUd1d2006h 



(2) AKARI; (3) IRAS; (4) This work; 



GAL and ATLAS GAL datfl and interferometric visibilities for 
comparison with the MIDI data. We do not attempt to match the 
near-IR VVV images, as AFGL 4176 is saturated in them. The 
relative goodness of fit of models was assessed by constructing 
a x 2 estimate which includes roughly equal contribution from 
these three measures. 



4.1 .1 . Constraints and uncertainties on fundamental 
parameters 

As noted in the introduction, the distance to AFGL 4176 cannot 
be reliably determined at present. We therefore chose to inves- 
tigate both near and far models, at distances of 3.5 and 5.3 kpc, 
respectively. 

Under the assumption that the H II region seen in radio ob- 
servations is powered by the AFGL 4176 source, the number of 
ionizing photons, and consequently a minimum inferred spectral 
type, can be dete rmined. The radio flu x measurements at 843 
and 1415 MHz o f lCaswell et al.l (Il992l) appear to be consistent 
with thermal bremsstrahlung radiation from an H II region. For a 
distance o f 5.3 kpc, using the m ain- sequence luminosities for O 
stars f romlMartins et al.l (120051) and the stellar atmosphere mod- 
els of K uruc zl (119791) . we find a minimum spectral type of the 
central source of about 08. This corresponds to an effective tem- 
perature of about 33 000 K, which we adopt for the temperature 
of the central blackbody source in our models. 

For the amount of foreground extinction, we adopt a simple 
prescription of 1 magnitude of visual extinction per kiloparsec, 
giving us Ay = 3.5 and 5.3 mag for the near and far models, 
respectively. This small amount of additional (i.e. interstellar) 
extinction was applied on top of our models, and has the effect 



For the PACS array on Herschel, we took 
the PSFs from the PACS calibration web site 
( http ://herschel. esac . esa.int/twiki/bin/view/Public/PacsCalibrationWeb ) . 
For the LAB OCA bolometer array used on the APEX telescope, we 
take t he PSF to be well- approximated by a Gaussian with FWHM = 
1972 (ISiringo et al.ll2009h . 



of further lowering the near-IR flux levels and slightly deepening 
the 10 yum silicate feature. 



4.1.2. Dust model 

One of the foremost practical constraints on the choice of dust 
model is the need for opacities at short wavelengths (A < 
0.1 /mi). Since the luminosity of this source is on the order of 
10 5 Lq and the amount of ionizing flux (inferred from the radio 
continuum observations) is equivalent to at least an 08 star star, 
the vast majority of the stellar radiation is emitted at UV and 
optical wavelengths. 

A number of attractive dust models for conditions expected 
in protostellar cores, envelopes a nd protoplanetary disks have 
been developed (for example by [O ssenkopf & Henningl 1 19941; 
iHenning & Stognienkol[l996l : ISemenov et al.1 120031) . However, 
a common shortcoming of these models is the lack of opac- 
ity information at short (UV) wavelengths, due to convergence 
problems of the grain models (D. Semenov, private communi- 
cation). On the other hand, the empirically derived opacities of 
lLaor & Drainel (1 19931) and subsequent works cover wavelengths 
from 100 A through 1 mm, with the caveats that the real physical 
makeup of the dust is unknown, and that at far-IR and (sub)mm 
wavelengths, where extinction is effectively impossible to mea- 
sure, the dust opacities are not well known. 

We therefore adopt a multi-component, temperature- 
dependent dust model, guided by the presence of known spec- 
tral features in the observations, and limited by the available 
coverage and reliability of different opacity models at different 
wavelength regimes. The transition between different regimes, 
and therefore the presence of different species throughout the 
surrounding material, are governed by the temperature structure 
iteratively solved by MoDust. 

Based on the presence of the 3.1 /mi water ice absorption 
feature in the ISO SWS spectrum, we infer that at least some 
of the amorphous silicate grains are surrounded by ice mantles, 
formed in the cold material of the parent molecular cloud. To 
include t his in our model, we use opacit ies from the coagulation 
model of lOssenkopf & Henningl (Il994l) . Following the sugges- 
tions in that paper, we interpolate between columns 1 and 2 of 
their Table 1 using the depth of the 3.1 jum ice feature. From the 
ISO SWS spectra, we measure T3.1 = 0.4, and use the resulting 
interpolated opacities for regions where T^ ust < 130 K. 



For grains with T^ ust > 13 K, we use the "astr onomical 
silicates" and graphite grains of lLaor & Dr aine (1993) with the 
"MRN" distribution of grain sizes (n(a) oc a~ 3 5 for grains be- 
tween 5 nm and 250 nm, where n is the number of grains, and a is 
the grain radius), which is generally appropriate for diffuse ISM 
conditions. We adopt a ratio of 60:40 for the silicate to graphite 
mass abundance in this temperature regime, although we note 
this is not directly constrained by observations. For the graphite 
grains, we adopt a destruction temperature of 950 K, which cor - 
responds to the process of OH sputtering (iDuschl et al.l 1 19961) . 
For the silicate grains, we do not enforce a destruction tempera- 
ture; at the inner edge of our envelope models, the silicate grains 
typically reach a temperature of about 2000 K. 

The dust model therefore consists of a region of "cold" 
(Tdust < 130 K) amorphous silicate and carbon grains with ice 
mantles, as well as "warm" (130 K< 7d us t < 950 K for graphites, 
and 130 K< 7d us t for silicates) amorphous grains similar to the 
"standard" diffuse ISM mix. 
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4.1 .3. Density structure 

For the density structure of the dust in our model, we adopt a 
piecewise power law consisting of concentric, spherical shells, 
with no gaps in between the shells. In order to simultaneously 
reproduce the observed visibility levels, far-IR/sub-mm spatial 
distributions and SED with this model, we find it necessary to 
include at least three shells, with a jump in density at the in- 
terface between the first (innermost) shell and the second shell. 
At the interface between the second and third shells, however, 
we were able to achieve acceptable fits without a density jump, 
thereby avoiding the need for an additional model parameter de- 
scribing this. We therefore characterize the density of material 
in the outer shells by the density p 2 3 at the contact point r 2 3 of 
the second and third shells. The precise structure used, therefore, 
has in principle nine free parameters (r\ 9 r 2 , r 2 ^, r , p\, P23, p\, 
P2, P3) and follows the form 



p(r) = 



Pi 



r\ < r <r 2 



P23( £ f) r 2<r<r 23 
P23( ! f) P3 r 23 <r<r 



(1) 



4.1.4. Results 

We present two fits to the observational data using the three- 
shell model: one for the "near" distance of 3.5 kpc, and one for 
the "far" distance of 5.3 kpc. The parameters are summarized in 
the first two lines of Tabled comparisons with the observations 
are shown in Figs. and [TO) the top panel shows the modeled 
flux levels from the photometric data from Table [2l the center 
panel shows the modeled MIDI visibilities at three wavelengths; 
and the lower panel shows the radial distributions from Herschel 
(70 and 160 /mi) and APEX (870 /mi). 

We defer a detailed discussion of the results to Sec. [5j 
but note the most prominent points of the model fits here. 
Both models reproduce the SED reasonably well at wavelengths 
longer than about 8 /mi, but significantly underestimate the near- 
infrared flux levels. Both models also reproduce the overall vis- 
ibility levels, as well as the spatial frequency of the zero cross- 
ing, however the far model provides a much better match, par- 
ticularly at the two shorter wavelengths (red and green points 
in the middle panels of Figures l9l and [TO. Finally, despite the 
good agreement of the radial profiles at wavelengths of 70 and 
870 /mi, both models overestimate the extent of emission at 
160 /mi, although the near model is in better agreement with the 
observed radial profile. The 160 /mi filter might be contaminated 
by strong [CII] emission at 157 /mi, detected in emission in the 
ISO LWS spectrum, which could alter the spatial profile of the 
brightness distribution seen in this broadband filter. However, 
since no spatially-resolved [CII] spectral map is available for 
AFGL 4176, this ad-hoc explanation is not beyond doubt. 

4.2. Two-dimensional geometric modeling 

While the one-dimensional approaches used above can provide 
useful insight, the deviations from spherical symmetry observed 
in the MIDI data should be explored. Such deviations can be- 
come especially clear in cases where multiple uv points have 
been measured with similar projected baselines (i.e., spatial fre- 
quencies), but differing position angles. Any significant varia- 
tions in the visibility amplitude for a given spatial frequency 
will thus require deviations from a circularly- symmetric inten- 
sity distribution to be explained. 
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Fig. 9. Model fits of the observational data for the near model. 
Top: The SED, where the aperture-corrected model flux levels 
are shown as crosses and the dotted line, and the observed flux 
values are shown as filled squares and the solid line. Middle: The 
modeled (lines) and observed (points) MIDI visibilities for three 
wavelengths. Bottom: The modeled (blue points) and observed 
(green points) radial intensity distributions from the Hi-GAL (70 
and 160 /mi) and ATLASGAL (870 fim) images. 



We see clear signs of such behavior for projected baselines 
longer than about 25 m, particularly for spatial frequencies larger 
than ~ 10 fringe cycles/arcsec. For small variations in spatial 
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Table 3. One-dimensional radiative transfer model parameters 



Model 


D 


r\ 


ri 


r 2 3 


r 


Pi 


P23 


Pi 


Pi 


P3 


M dust 


U 




(kpc) 


(AU) 


(AU) 


(AU) 


(AU) 


(g cm" 3 ) 


(g cm" 3 ) 








(Mo) 


(10 5 L ) 


Near 


3.5 


40 


65 


70000 


125 000 


1.71 x 10- 1 * 


3.56 x 10~ 22 


1.8 


1.05 


0.5 


4.70 


1.44 


Far 


5.3 


58 


108 


56000 


195 000 


9.22 x 10" 19 


4.71 x 10" 22 


1.2 


0.9 


0.5 


15.8 


2.88 


G91 


4.0 




100 


2646 


50134 




1.49 x 10" 21 




2.0 


0.0 


1.32 


1.57 



Table 4. Two-dimensional geometric fit parameters 



Parameter 


Best fit value 


r Q 


158 mas 1 


P 


0.467 


$ 


158° 


e 


59.7° 


F G IF D 


1.64 


g 


143 mas 2 



Notes. 1 552 AU at the near distance; 836 AU at the far distance. 
2 500 AU at the near distance; 757 AU at the far distance. 

frequency, we see visibility levels varying by over a factor of 5. 
Furthermore, there is a clear dependence on position angle: for 
measurements with V < 0.2, the visibility amplitude decreases 
with increasing position angle, while for measurements with V > 

0. 2, the visibility amplitude decreases with increasing position 
angle. 

This behavior (visibility amplitude increasing with position 
angle for V < 0.2, and decreasing with position angle for 
V > 0.2) cannot be modeled with a simple elongated Gaussian. 
However, such a trend is possible in the case of a source distribu- 
tion which has a well defined "edge" in its intensity distribution, 
causing a zero-crossing in the complex visibility with substantial 
power in the first side-lobe, for example a disk model. 

To model this asymmetric behavior, we use a simple param- 
eterized disk to match the dependence on position angle seen 
in our MIDI visibility measurements. We approximate the disk 
with a geometrically thin, optically thick blackbody emitter, with 
a temperature depending on radius in the form of a power law in 
the form of 

T (r) = T f (^) P n<r<r , (2) 

where Tf is the temperature at a fix-point ry, and r t and r are the 
inner and outer disk radii, respectively. We adopt Tf = 1500 K, 

1. e. the approximate evaporation temperature of silicate dust, and 

let rf - ^jL/4ncrTj, i.e. the radius at which a blackbody of 

luminosity L will have a temperature of 1500 K. For the lu- 
minosity of our far model, rf = 37 AU. We find that models 
in which the inner disk radius is at rf provide poor fits to our 
MIDI data, and a substantially smaller inner radius is required. 
We adopt r t = 1 AU, which is well below the spatial frequencies 
probed by our MIDI measurements. We thus implicitly assume 
the inner disk, consisting of gas, remains optically thick. This is 
in line with interferometric measurements of high-mass young 
stars, whose disk emission appears consistently more compact 
than expected for dusty, disk-only models (e.g. iMonnier et al.1 
l2005h . 

Finally, to account for the more symmetric behavior at 
shorter baselines, we include a one-dimensional component in 
the image intensity distribution in the form of a Gaussian cen- 
tered on the disk. Thus, as fit parameters, we use the outer radius 



of the disk r , the disk inclination 6 and position angle 0, the ex- 
ponent p of the temperature power-law, the flux ratio FqIFd of 
the Gaussian component to the disk component, and the FWHM 
6q of the Gaussian component. 

4.2.1. Results 

The parameters that best fit the observed MIDI visibilities were 
derived from a brute-force grid search of parameter space, where 
we further refine the parameters of the best fit from the grid 
search by using the downhill simplex algorithm. We find that this 
disk model produces a very good fit to our data, and in Table ffl 
we summarize the parameters of our best-fit model to the MIDI 
visibilities at a wavelength of 10.6 yum. We show the resulting 
model visibilities as triangles in Fig.[TTl and the model image in 

Fig.ua 

The observations are best described by a disk spanning up to 
several hundred AU, with the disk temperature decreasing grad- 
ually as a function of radius (T oc r -0 47 ). The disk is inclined 
(6 = 60°), and the semi-major axis is oriented 158° east of north. 
Besides the disk component, a substantial halo is also present, 
with a FWHM size equal to about half the diameter of the disk, 
and a total flux at A = 10.6 jim of 1.6 times the disk component. 



5. Discussion 

5.1. Comparison with previous models 

In order to gauge the importance of including spatial information 
in the modeling process, we reconstructed the one -dimensional 
radiative transfer model for this object presented by lGurtler et al.l 
(I1991I) . using their dust composition and density structure to- 
geth er with the M o Dust code used for our models. The model 
from lGurtleretaD (Il99ll) . which we refer to as G91, is a two- 
zone power-law model, and can be described using the density 
prescription from Eq.[T]if the parameters r\,p\ and p\ are omit- 
ted. In Tabled we show the parameters for the G91 model using 
this prescription. 

In Fig. [13] we show the resulting SED, mid-infrared visi- 
bilit ies, and far-infrared/ sub-mm radial intensity profiles. Note 
that lGiirtler etal.1 (Il99lf) fit only the SED to the extent it was 
known at the time, and that they did not have detailed spa- 
tial information (i.e. mid-IR interferometry and far-IR imag- 
ing) available. We indicate the measurements they had avail- 
able by squares in Fig. [T3J while the more recent Hi-GAL and 
ATL AS GAL measureme nts, as well as the 1.2 mm measurement 
from lBeltran et al.l (120061) . are indicated by diamonds. Regarding 
the flux values available at the time of their study, we note in 
particular that their 1.3 mm measurement of 131 mJy was prob- 
ably in error, as the position indicated in their Table 2 is more 
than 50" from the in frared source. More recently, for example, 
iBeltran et~aD (120061) reported a 1.2 mm flux value of 3.86 Jy 
for the source, which seems to be in good agreement with the 
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Fig. 12. Image of the fit of the disk model to the mid-infrared 
visibilities at A = 10.6 yum. 



Fig. 10. Model fits of the observational data for the far model. 
Figure notation is the same as for Fig. [9] 



ATLASGAL flux level of 1 1.56 Jy at 870 fim which we present 
here. 

In Fig.[l4l we show the density laws used for both our near 
and far models, as well as the G91 model. To account for the 
different distances to the source used in each model (see TableO, 
we express the dust density as a function of angular position on 
the sky. We also overlay the effective angular scales probed by 
our MIDI measurements, where have used the same color coding 
and grouping by projected baseline as in Fig.[T] 



For the near and far models, we see that the inner shell with 
enhanced density is probed primarily by the longest baseline 
measurements, while shorter baselines probe successively larger 
scales, out to about 071. Because this inner region is completely 
absent in the G91 model, the mid-IR visibilities at longer base- 
lines are vastly under-predicted by this model. This is despite 
the relatively good fit to the SED (with the exception of the sub- 
mm/mm regime), as well as the 70 and 160 yum radial profiles, 
which only emphasizes the importance of including spatial in- 
formation at different scales into such radiative transfer models. 
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Fig. 13. Fit to model by iGurtler et al.l (1199 lh . Figure notation is 
the same as for Fig. but in the top panel of the SED we dis- 
tinguish between flux values availab le today (filled diam onds) 
and those available at the time of the lGurtler et al.l (Il99ll) study 
(filled squares). 



5.2. The disk around AFGL 4 1 76 

The mid-infrared interferometric observations presented in this 
paper provide the first evidence for any signs of asymmetry in 
this deeply-embedded MYSO to date. Using a model consisting 
of a symmetric halo and an optically thick, geometrically thin 
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disk, we find that a disk with an inclination angle of 60° and a 
position angle of the semi-major axis of 158° east of north pro- 
vides an excellent fit to our MIDI data. The best-fit temperature 
profile of the disk is T oc r~ 0A1 . At longer wavelengths, the de- 
viation from spherical symmetry in our MIDI data is much less 
pronounced, which suggests that the emission from the disk is 
dominated by a more spherical halo component at these wave- 
lengths. 

Nevertheless, because of the wavelength range, mid-infrared 
observations probe primarily cooler material (~ 200 - 600 K), 
expected to be well away from the hotter, inner circumstel- 
lar regions where signs of accretion onto the central star may 
be present. Furthermore, the inner radius of the disk around 
AFGL 4176, adopted in this work to be 1 AU, cannot be de- 
termined by our observations, as it remains unresolved by the 
MIDI data. In this context, near-infrared interferometric obser- 
vations which probe these inner regions would make a superb 
compliment to the current data set. 

On the other hand, the lack of detection of an outflow com- 
ponent, or any other larger-scale asymmetries, remains an open 
question. This may, for example, suggest that the object is still at 
a very early stage in the formation process, and that large-scale 
outflows have not had time to form. In any case, far-IR/sub-mm 
interferometric observations with ALMA may provide the key to 
bridging the symmetric structure on scales of tens to hundreds of 
thousands of AU to the asymmetries observed at tens to hundreds 
of AU. 



5.3. Comparison with other MYSOs 

To date, very few MYSOs have been studied on scales of tens 
to hundreds of AUs. Here we compare our results with similar 
studies of five other MYSOs. In Table[5]we list the luminosities, 
the depth of the 10 yum silicate absorption feature, whether or not 
any compact H II region is present, and whether or not any out- 
flows have been detected, which can be useful for independently 
constraining disk geometry. As previously noted, however, in the 
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Table 5. Summary of parameters for several MYSOs studied on scales of ~ 10 - 100 AU. All values approximate. 



MYSO 


D (kpc) 


L(L Q ) 


A v (mag) H II region present? 


Outflows present? 


References 


AFGL 4176 


3.5-5.3 


10 s 


100 


Yes 


? 


1,2,3 


W33A 


3.8 


10 5 


280 


Yes 


Yes 


4 


NGC 3603 IRS 9A 


7 


10 5 


? 


Yes 


? 


5 


AFGL 2136 


1.7 


2 - 3 x 10 4 


100 


Yes 


Yes 


6,7 


IRAS 13481-6124 


3.2 


3xl0 4 


40? 


Yes 


Yes 


8 


M8E-IR 


1.5 


7x 10 3 


30 


No 


Yes 


9 



Refer ence s. (1) This wor k; ( 2)lPhillips et all (fl998h : (3) lEllingsen et all (120051) : (4) Ide Wit et all (12010) : (5) IVehoff et all (12010) : (6) Ide Wit et all 
J20TT1) : (7) lFollertl(l20TTI) : (8) lKraus et all d201Qh : (9) lLinz et all (12009^ 7 



case of AFGL 4176, no outflow compon ents have been detected 
dDeBuizer 20031: iDe Buizer et al 1 12009b . 

The sources W33A and NGC 3603 IRS 9A have similar lu- 
minosities to AFGL 4176, on the order of 10 5 L©, and both have 
associated hyper-compact H II regions. However, IRS 9 A seems 
to be much less embedded: besides having an optical component, 
the Af-band spectrum shows no 10 /mi silicate feature at all. At 
mid-infrared wavelengths, the source is very extended, and com- 
pletely overresolved in M IDI observations probing ~ 200 AU 
scales (IVehoff etal.1 l2010h . W33A, on the other hand, could 
be even more em bedded than AFGL 4176 (Ay ^ 200 mag, 
Ide Wit et all 1201 Ol) . but the Af-band visibility levels are poorly 
determined due to the low amount of (correlated) flux in the 
de ep 10 £mi silicate abs orption feature. The source was modeled 
by Ide Wit et all (120101) . who concluded that most of the mid- 
infrared emission on scales of 10 2 AU is dominated by colli- 
mated outflow cones, without compelling evidence for signifi- 
cant emission from a dusty disk. 

At an estimated luminosity of ~ 5 x 10 4 L©, the SED 
of AFGL 2136 closely resembles that of AFGL 4 176, includ- 
ing a similar depth of the 10 fim silicate feature (ISmith et al.l 
2000 |). Usinff visibi lity measurements obtained with MIDI, both 



de Wit et al.f(l201ll) and lFollertl(l201l[) repo rted mid-IR emission 
scales of ~ 200 AU, while IFollertl ( 201 lb found an elongated 



structure at these scales, parallel to earlier predictions of disk 
orientation. 

Despite having a lower luminosity, on the order of 3 x 
10 4 L Q , the SED of IRAS 1348 1-6 124 is somewhat similar to 
that of AFGL 4176, although th e 10 /im silic a te fea ture is much 
shallower in IRAS 1348 1-6 124. iKraus et al.l (l2010h discovered 
a molecular outflow, and used near-IR interferometric obser- 
vations to resolve an elongated structure oriented perpendicu- 
lar to the outflow direction. The extended outflow shocks de- 
tected at near-IR wavelengths at large distances from the source 
(Stecklum et al. 201(3), as well a s the radiative transfer mod- 
els presented by lKraus et al.l (120101) . suggest the disk is strongly 
inclined. However, the 20 M^ (gas + dust m ass) disk at scales 
of 5-150 AU proposed by lKraus et al.l (120101) is clearly a funda- 
mentally different model than the structures we present here for 
AFGL 4176. For comparison, the dust mass enclosed in the near 
and far models for r < 150 AU is 8 x 10" 6 M© and 4 x 10" 5 M©, 
respectively, although a direct comparison between studies using 
near- and mid-IR interferometric measurements is difficult. 

Finally, the source with the lowest luminosity, M8E-IR, has 
an overall extinction of Ay « 30 m ag, as determined from the 
depth of the 10 //m silicate feature dLinz et al .1120091). Outflow 
signa tures have been detected for this source dMitchell et~aD 
Il988l) . No cm continuum emission has been detected directly 



coinciding with the IR object. This lack of an ultra-compact H II 
region might be related to th e hypothesis of a bloated, cooler 
central object, as suggested by lLinz et al.l (l2009h . 

Given the small selection of sources for which this infor- 
mation is available, it is not possible to infer clear, empirical 
correlations between the properties listed in Table \5\ Of the 
sources presented, none of them seem to be direct analogues 
of AFGL 4176, however W33A might be the closest match in 
terms of luminosity, SED a nd emission scales at m id-infrared 
wavele ngths. Suggestions bv lde Wit et al] (1201 lh and lLinz et al.1 
rf20T!h indicate that mid-IR interferometric observations of at 
least some MYSOs may be heavily influenced by outflow com- 
ponents. On t he other hand, the mid-IR visibilities presented by 
IFollertl (120111) apparently trace the disk structure of AFGL 2136. 
Taken together, it becomes clear that the mid-IR emission struc- 
ture of MYSOs at scales of tens to hundreds of AUs can be com- 
plex, and caution should be exercised in interpreting data sets 
with limited uv coverage. 

6. Summary and conclusions 

In this paper, we reported on spatially resolved observations 
and modeling of the massive young stellar object AFGL 4176. 
We presented the largest amount of Af-band visibilities for any 
MYSO to date, obtained with MIDI on the ESO/VLTI, together 
with maps of extended emission at 70 and 160 /jm from the Hi- 
GAL survey on the Herschel Space Observatory and 870 /im 
imaging from the ATLAS GAL survey with the APEX telescope. 

With the exception of the mid-infrared interferometric mea- 
surements presented in this paper, the observational data avail- 
able for this deeply-embedded object do not provide indications 
of asymmetry or preferred geometry. In view of this, we used 
spherically symmetric radiative transfer models, consisting of an 
envelope with a density enhancement at its inner edge, to model 
the source. Using these one-dimensional models, we were able 
to satisfactorily reproduce, simultaneously, the SED and spatial 
structure at mid-IR through mm wavelengths. 

However, the one-dimensional models by definition cannot 
reproduce the asymmetries present in the MIDI data. We there- 
fore modeled these data separately, using a multiple-component 
geometric analysis. We interpret the MIDI observations in terms 
of a circumstellar disk around AFGL 4176. If this disk hypoth- 
esis is correct, it would make AFGL 4176 the most luminous 
(massive) young star to have a disk detected around it to date. 
However, beyond initial detection, much remains to be clarified 
regarding the nature and extent of the disk, making this object a 
top priority for additional interferometric imaging, at both near- 
IR and far-IR/sub-mm wavelengths. 



http://www.jcu.edu.au/hmsfl0/Presentations/Monday/Session3/Stecldum.pdf 
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